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ABSTRACT 

We present a search for long-period variable (LPV) stars among giant branch stars in 
M15 which, at [Fe/H] ~ -2.3, is one of the most metal-poor Galactic globular clusters. 
We use multi-colour optical photometry from the 0.6-m Keele Thornton and 2-m 
Liverpool Telescopes. Variability of SV ^ 0.15 mag is detected in K757 and K825 over 
unusually-long timescales of nearly a year, making them the most metal-poor LPVs 
found in a Galactic globular cluster. K825 is placed on the long secondary period 
sequence, identified for metal-rich LPVs, though no primary period is detectable. 
We discuss this variability in the context of dust production and stellar evolution at 
low metallicity, using additional spectra from the 6.5-m Magellan (Las Campanas) 
telescope. A lack of dust production, despite the presence of gaseous mass loss raises 
questions about the production of dust and the intra-cluster medium of this cluster. 

Key words: stars: AGB and post-AGB — stars: late-type — stars: Population II — 
stars: variables: other — stars: winds, outflows — globular clusters: individual: M15 



1 INTRODUCTION 



The onset of radial pulsation on the red giant branches is 
important in stellar evolution: it is one of several linked phe- 
nomena (including dust production and substantial mass 
I loss) that control the endpoint of stellar evolution and injec- 
, tion of mass into the interstellar medium. While these phe- 
nomena are linked, the relative timing of their ap pearance is 
poorly understood (|McDonald fc van Loorj[2007^ . In partic- 
ular, it is debatable whether pulsation can provide enough 
energ y to assist mass loss in these evolved stars (e.g. iBowenI 
Il98; 



Optical photometric variability in highly-evolved stars 
is known to be less p ronounced in metal-poor systems. 
iFrogel fc Whitelockl (fl998 ) considered Long-Period Variable 
stars (LPVs) in both the Galactic Disc and globular clus- 
ters in order to prove this dependence, though they only 
consider stars showing large-amplitude variability as LPVs 
and do not consider semi-regular variables (SRVs, which are 
included in the definition of LPVs for the purposes of this 
work). This raises the question of whether pulsation is ca- 
pable, or required, to drive mass loss from metal-poor stars. 

We describe here the search for LPVs in one of the most 
metal-poor Galactic globular clusters, M15. This cluster is 
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the only one known to harbour a dusty and/or gaseous in- 
terste llar medium and has several infrared-excessive giant 
stars (|Evans et al.ll2003l : Ivan Loon et al]|2006bl : iBover et al.l 
l2006l ). giving the strong implication that dusty stellar winds 
are present in the cluster and, by further implication, 
pulsation-driven winds. 

No LPVs have so-far be en found in M15, though several 
candidates were identified bv lMoslev fc White! l| 19751 ) (Table 
[TJ. These were followed- up by I Welt^^ ( 19851 ) . who could not 
find photometric variability in three targets (K169, K288 
arrd K709) and retain ed K757 and K 825 as candidate vari- 
ables (identifiers from lKustner|[l92ll ) . 



2 OBSERVATIONS 
2.1 Liverpool Telescope 

The first dataset for thi s work comes from the two-metre 
Liverpool Telescope (LT: ISteele et al.ll2004l ). situated on La 
Palma. A total of 38 observation blocks were taken, span- 
ning 462 days, from 2007 April 23 to 2008 July 27. A gap is 
present between 2007 December 04 and 2008 April 23 when 
the cluster was viewed in too close proximity to the Sun. Ob- 
servations were taken roughly every eight days, with the ex- 
ception of a one-week block of daily observations in August 
2007, to build in redundancy against shorter-period variabil- 
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ity. Scheduling and weather constraints mean the observa- 
tions are randomly distributed enough not to form strong 
aliases in a Fourier spectrum. 

Each observation consisted of 3 x 6 s 3' (477 nm) and 
3 X 5 s i' (762.5 nm) exposures, thus minimising stochastic 
effects from cosmic rays, etc. An additional 2 s image was 
taken in both filters in case the stars saturated the detec- 
tor. Preliminary analysis showed this was not necessary and 
these were not taken in the 2008 observing season. Image 
calibration was performed automatically before receipt. 
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2.2 Keele Thornton Telescope 

Additional images were taken with Keele's 24-inch (60-cm) 
Thornton Telescope (KT) sited on Keele University campus, 
with elevated (205m above sea level) views over the Stafford- 
shire countryside. Although these images have poorer pho- 
tometric accuracy and astrometric resolution, they cover a 
longer timebase. The data were taken using the Santa Bar- 
bara Instrument Group ST7 CCD camera (at Newtonian 
focus). Several features were found that have the potential 
to affect the photometry: namely tracking, coolant issues 
and random bias and dark current problems. In practice, 
however, we find that the statistical scatter in photometry 
from this telescope is smaller than that from the LT. 

With these issues in mind, observation blocks were de- 
vised of multiple 10-second observations, taken in V-, R- and 
/-bands (550, 700 and 880 nm, respectively). Observations 
taken at low altitude have not been corrected for differential 
reddening (due to its night-to-night variation): some of the 
scatter in the resultant photometry may be attributable to 
it. 

In total, 76 observation blocks were taken, spanning 
from 2006 January 14 to 2008 November 24, though most 
observations were taken after 2007 August 02. Weather con- 
straints at Keele are significantly worse than those on La 
Palma, leading to even less-regular sampling. 
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Figure 1. Negative example images of M15 from the Keele 24- 
inch Thornton Telescope (top panel; V-band image) and the 
Liverpool Telescope (bottom panel; average of g' and i' im- 
ages). Comparison stars (blue circle s) and candidate variables 
(red squares) from Mosley & Whit(| l ll975l ) are shown. The in- 
set in the Keele figure shows the heavy crowding in the core. The 
Liverpool Telescope figure is displayed with a gamma-law scaling 
to aid viewing. 



3 RESULTS 



3.1 Candidate selection 

The presence of variability can be determined via differ- 
ence imaging. Here, an image {A) is scaled and convolved to 
match the point-spread function of a second image [B, taken 
at a different epoch) . The resulting image is subtracted from 
image B to create a difference image. Variable stars show up 
as non-zero flux. 

We performed difference imaging on pairs of images 
from both telescopes, taken at various epochs. Variability 
was immediately identiflable in K825 and suggested in K757. 
Of the other stars in the cluster, the only obvious suggestion 
of variability came from the cluster's core, which is entirely 
unresolved in our images. We therefore do not consider vari- 
ab ility in stars other than those candidate variables listed 
in lMoslev fc White! l| 19751 ). which include K757 and K825. 



3.2 Photometric reduction 

3.2.1 Comparison star selection 

Example images from both telescopes are shown in Fig. [1] 
Aperture photometry was extracted from these ima ges us - 
ing the software AIP4WlN v. 1.3.5 (p3errv fc Burnelj 120051 '). 
Comparison stars (Table (2] Fig. [l| were chosen to be: 

• isolated (resolved from any visible neighbours by at least 
2.5 X the full- width half-maximum seeing), so that seeing has 
minimal effect on photometry due to blended stellar images; 

• bright (i' < 14 mag), to provide good signal-to-noise, but 
not saturated; 

• red (B — V > 0.75 mag), to avoid problems with differen- 
tial reddening at high airmass; 

• invariable, nor variable candidates themselves (we retain 
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Table 2. Comparison stars used in the photometric reduction of 
images from the Keele Thornton and Liverpool Telescopes. K497, 
K508 and K589 were not used in the final analysis due to possible 
intrinsic variability. 



Designator 



Co-ordinates (J2000) V {B - V) 



(1) 


(2) 


(3) 


(4) 


RA 


Dec 


(5) 


(5) 








Comparison stars (Keele) 






K87 


III 


438 


S29 


2129 45.81 


12 08 45.1 


13.81 


1.01 


K112 


II 


213 




2129 47.90 


121131.3 


15.12 


0.883 


K114 


III 


437 




2129 47.86 


12 08 45.0 


13.87 


1.06 


K129 


II 


187 




2129 48.61 


121145.6 


14.27 


0.96 


K144 


11-75 


248 




2129 49.78 


121105.6 


13.00 


1.25 


K158 


III-33 


414 




2129 50.27 


12 09 02.5 


14.13 


0.78 


K731 


1-63 






2130 00.52 


121136.9 


14.24 


0.96 


K912 


IV-48 


457 




2130 04.20 


12 08 27.4 


14.40 


0.91 


K1040 I 


319 


S6 


213010.49 


121006.2 


13.40 


1.18 



Comparison stars (Liverpool Telescope) 



K224 


n-30 


2129 52.30 


121051.3 


14.51 


1.07^ 


K238 


11-29 


2129 52.63 


121043.8 


14.46 


1.17^ 


K319 


II 


2129 54.42 


121102.4 






K497 


III 


2129 57.17 


12 0917.3 


14.42 


1.29^ 


K508 


III 


2129 57.49 


12 09 06.1 






K589 


IV 


2129 58.71 


12 08 55.8 






K928 


1-43 


2130 04.61 


121032.6 


13.83 


1.03 


K943 


1-38 


2130 05.14 


1210 04.1 


14.24 


0.94 


K961 


1-41 


2130 05.69 


121015.7 


14.09 


1.03 



^Kustneij l|l92lh : - ^Ard {TgH) (incomplete designations de- 
note the quadrant of Arp's map containing the star if the star 
is without a designation); ' ^Buonanno et aL I lll983h : ISandagel 
(|l970h : ICudworthI l|l976l); '^f rom photographic and photovi- 
sual magnitudes from iBrowrJ lll95ll ). 



K1040 as its known short-p eriod variability of 0.04 mag i s 
well below our sensitivity — iBao- Anlll990l : lYao et al.|[l99a ) . 



provided by AIP4Win) were calculated manually by adding 
in quadrature the relative Poissonian errors in the photo- 
metric signal and calibration frames. 

The absence of /ong-period variability of the compari- 
son stars K144 (Arp 11-75) and K1040 is interesting in itself, 
as these stars are identified as be ing among the brig htest, 
reddest variables in the cluster bv lFrogel et al.l (|l983l ). (Al- 
though K1040 (S6) is a known short-period variable ( i]3.2.H) . 
the 4.3-hour period is too short and the 5V = 0.04 mag am- 
plitude too small to affect our data.) 

An identical photometric reduction was done on the 
candidate variable stars. The flux of each of the candidate 
variables was divided by the sum of the fluxes of the com- 
parison stars to yield a differential magnitude. 



3.2.3 Liverpool Telescope 

A similar reduction was performed on the LT data. The fleld 
here is much smaller, though the seeing is generally better. 
A different set of comparison stars (Table [5} was used, fol- 
lowing similar selection criteria. This builds in redundancy 
against comparison star variation as we have two completely 
independent data sets. 

Fluxes and differential magnitudes were calculated 
on each processed image. These were summed over each 
night to improve fidelity. Here again, the non- variability of 
K238 (Arp 11-29) is interesting, as it is suggested to be a 
comparatively-bright AGB star (iFrogel et al.lll983l ). 

Internal comparisons between comparison stars show 
possible small-amplitude variations in K497, K508 and K589 
on the order of 20-30 mmag. The residual long-term varia- 
tion in the other six stars is <10 mmag. When constructing 
a summed comparison flux for each image, these three stars 
were discounted. 

The combined results from both the KT and LT are 
shown in Fig. [5] Differential magnitudes of the target stars 
are listed in Tables [3] &|1 



3.2.2 Keele Thornton Telescope 

The obvious problem is that bright, red stars tend to be 
other cluster giants, which may also vary. Our comparison 
stars show no intrinsic variation, though a systematic offset 
in several of the comparison stars (K129, K144 and K912) 
appears depending on the orientation of the telescope. The 
reason for this is unclear, though seeing and airmass may 
play a role: the telescope must be inverted at this declina- 
tion to observe at higher airmasses in the east, where hght 
pollution is greatest. There may also be a weak correlation 
between calculated differential magnitude and seeing qual- 
ity. Due to concerns that this offset is also present in the 
candidate stars, we have discounted the frames taken when 
the telescope was inverted from further analysis. 

Images were calibrated and stacked using AIP4WiN to 
increase signal-to-noise. Though not a lossless procedure, 
tests on individual and stacked images shows that the likely 
effect on the results is <1 milli-mag (mmag). A standard 
photometric reduction was performed: i.e. bias and dark 
frame subtraction, and flatfield division. Aperture photome- 
try was performed using AIP4WlN. Photometric errors (not 



3.3 Colour— magnitude diagrams 

Colour-ma gnitude diagr ams (CMDs) were generated using 
DAOPhot dstetsonll 19871 ). In contrast to the earlier aperture 
photometry, PSF-fitting was used to extract stellar magni- 
tudes from the images here in order to increase coverage 
toward the very crowded cluster centre. 

The errors calculated from PSF-fitting are correspond- 
ingly larger than those from aperture photometry: in the LT 
data, these are typically ±9-15 mmag for PSF-fitting, com- 
pared to ±1-2 mmag for aperture photometry. However, the 
aperture photometry errors do not take into account close 
stellar blends: the scatter of magnitudes for non-variable 
stars suggests that the actual errors in our aperture photom- 
etry may be larger, though still smaller than the PSF-fitting 
errors. 

After PSF-fltting, an offset was applied to the pixel co- 
ordinates from the i' image to match those in the g' image, 
and the closest object to each detection matched between 
the bands, within a maximum distance of flve pixels (1.3"). 
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Table 1. Candi date variab l e gian t stars in M15 from lMoslev fc White! l|l975l l with variability results. Proper motion memberships and 
colours are from lCudwortiil l ll976l ). Other references and notes as Table |2] Note that all periods are very approximate. 



Designator 




Co-ordinates (J2000) 


Membership v V 


{B- 


V) 


Variability 


Period 


(1) (2) 


(3) (4) 


RA Dec 


probability (km s~^) (mag) 


(ma 


g) 


found? 


(days) 



K134 


III-8 


416 


2129 49.07 


12 09 02.0 


~0% 




14.02 


0.99 


No 




K147 


III-34 


442 


2129 50.01 


12 08 43.4 


~0% 


+19-"^ 


12.62 


0.75 


Maybe 




K169 


11-64 


212 


2129 50.81 


121130.0 


60% 




13.48 


1.12 


No 




K288 


11-16 




2129 53.78 


1210 20.1 


90% 


-107.2^^ 


13.59 


1.06 


No 




K709 


IV-58 




2130 00.39 


12 07 36.0 


98% 


-100.7^ 


13.52 


1.07 


Maybe 




K757 


IV-38 




2130 00.91 


12 08 56.8 


95% 


-113.5*^ 


12.58 


1.36 


Yes 


~250? 


K825 


1-12 


S4 


2130 02.23 


121121.5 


67% 


-98.87* 


12.69 


1.37 


Yes 


350±20 



teustned 1I1921I '): ^Ard l ll955t) : ^Buonanno et al.l lll983l '): ' ISandagel lll970h : heliocentric radial velocities from ^Smoker et al.l 
1 I2OO2I) : iGebhardt et al.l lll997l '): iMeszaros et al.l ll2008al) : and ^Pilachowski et al.l llioool) ; cf. M15's velocity of tJhdioccntric = 
— 107.9 km s~^ and internal velocity dispersion of ~10 km s" -1 llGebhardt et al.lllQQTl V 
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Figure 2. Light curves of M15 candidate variables. Data have been offset vertically for clarity. Blue asterisks and magenta hollow squares 
represent Liverpool Telescope g' and i' photometry, while green crosses, orange filled squares and red dashes show V-, R- and 7-band 
Keele Thornton Telescope data, as labelled in the light curve of K134. 
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Table 3. Difference photometry of comparison stars with variable stars in data from Keele Thornton Telescope. A full electronic table 
is available online. 



Julian Apparent photometric magnitude 



Date 


K169 


K825 


K288 


K147 


K757 


K288 


K709 


V 


2453962.43 
3986.36 
3994.32 


1.408 
1.441 
1.462 


1.231 
1.223 
1.208 


2.385 
2.306 
2.321 


1.145 
1.109 
1.110 


1.464 
1.457 
1.449 


2.554 
2.722 
2.626 


1.211 
1.206 
1.182 



Table 4. Difference photometry of comparison stars with variable stars in data from Liverpool Telescope. A full electronic table is 
available online. 



Julian Apparent photometric magnitude 

Date K169 K825 K288 K147 K757 K288 



g' 



2454213.71 2.300 1.611 2.496 1.108 1.694 2.692 
4220.69 2.263 1.613 2.477 1.092 1.717 2.666 
4237.64 2.278 1.620 2.467 1.078 1.683 2.661 



Due to the lack of literature observations with g'-, i'- 
and i?-band magnitudes, and lack of good colour transfor- 
mations to obtain these from B- V- and /-band data alone, 
we have no absolute flux measurements to calibrate the 
(?'-, i'- and 7?- bands. To circumvent this, we have computed 
theoretical AB magnitudes of K238 in these bands, using 
the Sloan DSS filter transmission^ and the spectral energy 
distribution (SE P ) par ameter-estimation technique used in 
iMcDonald et al.l (|2009l ). to give an approximate physical 
magnitude. We find the star has a temperature of 4583 K 
and luminosi ty of 805 L©, assuming E{B ~ V) = 0.10 mag 
l|Harris|[l99^ '). and calculate g' = 13.783 mag, i' = 12.506 
mag and R = 12.532 mag. Note these values may be in error 
by -0.1 mag. The final CMD is shown in Fig. E] 

The CMD clearly resolves the red giant branch as a well- 
denoted sequence extending from a clearly-defined RGB-tip 
at i' ~ 12 mag, down to the main-sequence turnoff at i' « 19 
mag. The horizontal branch (mainly around i' ~ 16 mag) 
can also be seen extending to very blue colours. The dip at 
around (<?' — i') = —0.3 m ag has also been found in earlier 
studies (e.g. IStetsonlll994l ). as has the kink and gap at the 
very top of the giant branch {i' ~ 12.3 mag). Much of the 
scatter away from the giant and horizontal branches comes 
from poor magnitude determination or cross-matching near 
the crowded cluster centre. 



3.4 Variability 

Although variability is clearly seen in Fig.[2]in several stars, 
the strength, periodicity and regularity of variation is diffi- 
cult to determine due to the long-period nature of the vari- 
ations. 

In order to determine where variability occurs, we have 
^ http://www.sdss.Org/dr3/instruments/imager/#filters 
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Figure 3. Colour-magnitude diagram of M15 from the Liverpool 
Telescope for the night of 2007 October 03. The asymptotic and 
red giant branches (AGB/RGB) and horizontal branch (HB) are 
identified, the main sequence turnoff is at i' ^ 19 mag. 



created the index /i. We calculate this as follows: for each 
star and filter, we take the stellar magnitude time series and 
remove any points that deviate from the last observation by 
more than 0.15 mag. This removes any bad data deviating 
by 5J 5x the standard photometric error from the average 
observation. For this clipped dataset, we take the standard 
deviation, as and divide it by the standard deviation of the 
point-to-point differences, ap, giving /x = as/ap. For a non- 
variable star, n should be very near-unity. Slight departures 
from unity may be due to remaining bad data, but values 
significantly above unity should be due to variability. The 
index is equally sensitive to both regular and irregular vari- 
ability. 

For each star, we have listed in Table [5] two variability 
statistics: for the KT, given by the average of /i in the 
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Table 5. Variability indices of candidate variables, for the KT 
(^k) and LT respectively. Values of > 1 suggest vari- 

ability, the derivation is described in the text, "r.m.s." denotes 
the average root mean square variation over all filters from the 
KT (V, R, I) and LT (g' , i'), respectively. The value r is the 
correlation co-efficient between the KT and LT data. 



ID 






r.m.s. 

(mmag) 


r 


Variability 
detected? 


K134 


1.19 




30 






No 


K147 


1.36 


2.23 


15 


24 


-0.35 


Indefinite* 


K169 


1.14 


0.94 


18 


20 


-0.36 


No 


K288 


1.28 


1.33 


31 


11 


-0.40 


No 


K709 


1.35 




19 






Maybe 


K757 


1.65 


2.13 


44 


46 


+0.73 


Yes 


K825 


1.51 


1.98 


34 


41 


+0.91 


Yes 



* See ^4T2l for details. 



V-, R- and /—bands; and for the LT, which is the average 
of IJ, in the g' and i' bands. We also give the r.m.s. varia- 
tion, which demonstrates the stability of the KT data: the 
only star with a higher r.m.s. in these data is K288, which 
is close to the core and suffers from some blending. Finally, 
we also include the Pearson product-moment correlation co- 
efficient, r, between binned KT and LT data (positive values 
approaching unity denote good correlation). Here, each tele- 
scope's data is represented by a 30-day bin of the deviation 
from average magnitude, over all filters. This assists differ- 
entiation between intrinsic variability and variability caused 
by changing observing conditions. 



From Table O we can see that both telescopes show 
high values of — indicating the presence of variability — 
in K757 and K825. Marginal detections of variability are 
possible in K709, for which there is no LT data, and K147. 
In the latter case, however, there appears no correlation be- 
tween the trends observed from the KT and those from the 
LT, so we neither claim nor refute variability here. 



4 DISCUSSION 
4.1 Individual stars 

The results for individual stars are summarised in Table [T] 

4.1.1 K134 (Arp III-8) 

Although not always in the field of view of the LT, there is 
nevertheless sufficient data from both telescopes to make two 
independent tests of variability in this star. There appears 
to be no variation above the scatter in the data, and we can 
thus conclude that K134 does not show any substantial long 
period variations . Note that this star is not thought to be a 
cluster member (|Cudworthlll97^ '). 

4.1.2 KI47 (Arp ni-34) 

One of the brightest stars in the c luster, K147 is a relatively 
blue star: iBuonanno et al.l l|l983l l show it to have a near- 



zero {U — B) colour. Its proper motion and radial velocity 
show it to be a field star (Table [l|. While our observations 
show some variability may be present ( H3.4p . data from the 
two telescopes are discordant about the trend of this vari- 
ability, suggesting the observed variability is spurious. It is 
possible that these variations in our data are due to seasonal 
instrumental effects. The sensitivity of our data means that 
we cannot rule out variability in this star at low (J; ±1%) 
amplitudes. 

4.1.3 K169 (Arp 11-64) 

The photometry of K169 has a relatively low scatter. It 
clearly shows no large-amplitude variations. Our photome- 
try does not suggest variability in this star and we therefore 
conclude that it does not pulsate significantly on timescales 
from days to up to a year. 

4.1.4 K288 (Arp 11-16) 

K288 shows relatively low-noise photometry in the LT data, 
but shows a larger scatter in the KT p h otome try due to 
its proximity to the cluster core. IWeltvl (|1985| ) could not 
determine variability in this star because of this, though 
our data show low enough scatter to determine that there 
appears to be no coherent variation in this star of 5V ^ 0.01 
mag. We conclude that this star is also not variable on long 
timescales. 

4.1.5 K709 (Arp IV-58) 

K709 may show variability of 5V ^ 0.08 — 0.10 magnitudes, 
with maxima at JD 2 454 400 and 2 454 600, and a minimum 
near JD 2 454 500. The star also has a moderately-high vari- 
ability index in the KT data, comparable with that of K147 
( ^3.4|l . We do not have data from the LT to determine if 
this variability is real (it lies beyond the LT field of view), 
and the variations are close to the noise level in t he KT 
data. Variability in this star would be interesting, as IWeltvl 
(1985) find this star to be non-variable, with an r.m.s. of 
0.054 mag: the variations seen in our data would seem to be 
slightly larger. It is possible that Welty missed variability 
here due to sparsity of data coverage. 

Fig.|3]shows a Fast Fourier Transform (FFT) in each fil- 
ter. The left- most peak near cycles day"'^, most obvious in 
K709 in the R-h&nA, is the timebase of our photometry, but 
K709's power spectrum shows no other peaks of statistical 
significance. The data therefore appear to be of insufficient 
coverage to determine any periodicity, should it exist. 

4.1.6 K757 (Arp IV-38) 

As one of the brightest and reddest stars in the cluster, 
and also as a confirmed radial velocity member (Table [1]), 
this star is a prime candidate for being an LPV. Indeed, 
the data from both telescopes shows a concurrent dimming 
around JD 2 454 340 for ~100 days. This drop is wavelength- 
dependent, showing a clear increase with decreasing wave- 
length, decreasing from ~0.2 magnitudes for the g' data to 
<0.05 magnitudes in J-band. Further, smaller drops may be 
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seen near JD 2 454100 and JD 2 454 600, though they are 
less-well covered. 

The light curve is very similar to a typical binary star 
showing a single eclipse and ellipsoidal variations, though 
the star would have to be in a contact or near-contact binary 
to show variation at this (approximate) period. We consider 
it more likely that this star is showing pulsation. 

Again, the data are of insufficient coverage to deter- 
mine whether a regular period exists to this variability. The 
FFT (Fig.|4]) shows a broad peak in all bands (except i') be- 
tween 100 and 500 days, though visual examination of the 
lightcurve (Fig. ^ suggests that it is probably ~250 days. 
Further data are needed to see if this variation is truly pe- 
riodic. 



4.1.7 K825 (Arp 1-12) 

K825 is another particularly bright and red star in the clus- 
ter which is also a confirmed radial velocity member (Table 
[l]). Physically, the star is also the most-luminous star in 
the cluster, and one of the coolest. While we cannot dis- 
tinguish between AGB and RGB stars at this luminosity, 
K825 is likely to define t he RGB tip. It is also a known Ha 
variable (see ^4.2.21 and iMeszaros et al1l2008bh . suggesting 
either chromospheric activity or pulsation. 

Variations from this star were immediately visible in 
the data from an early stage of the analysis. A clear dip 
is seen near JD 2 454 300 and another near JD 2 454 650 
in all five bands. A further possible dip at JD 2 454 050 
may be present, but there is only one epoch of observation 
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here. Again, these dips are wavelength-dependent, decreas- 
ing from ~0.17 mag in g' to ~0.08 mag in 7-band. 

The hght curve shows a classic long-period variation at- 
tributable to stellar pulsation, with a rapid drop in bright- 
ness at bluer wavelengths as the star expands and cools, 
then a slower recovery as the star contracts and warms. FFT 
analysis (Fig. [Jjl shows peaks in all bands between 330 and 
370 days. The strong peak in the LT data at 200 days is an 
alias of this caused by insufficient data coverage. The lack 
of a clear drop at JD 2 453 950 suggests that the amplitude 
of variability is irregular and thus that this star is an SRV. 



4.2 Variability in the context of stellar evolution 

4-2.1 Pulsation as the source of variability 

We have found variability in two of the cluster's most 
evolved members (K757 and K825). These stars seem to 
share the metallicity of the cluster: both are within 0.12 
dex of the the clu ster's average value, [Fe/H] = -2.4 
l|Sneden et al.l Il997l . 2000) . These therefore represent the 
most metal-poor long-period variables known in our Galaxy. 
Outside our Galaxy, there is only an unconfirmed member 
of the Bootes I dwarf galaxy that may be more metal-poor 
llDairOra et al.ll2006l ). This star has a period of P 85 
days, SV « 0.4 and an optical magnitude which would place 
it on the AGB, w ere it con fi rmed a member. 

Notably, Sne den et al.l (|2000l ') also find both K757 and 
K825 to be barium-enhanced by [Ba/Fe] ~ -fO.37, among 
the highest enhancement in the cluster, su ggesting convec- 
tive m ixing of s-process elements (see also van Loon et al.l 
|[2003))- In Sneden et al.'s earlier work, they further suggest 
that both may be double stars, on the basis of "weak, blue- 
displaced extra abs orption components". T wo-epoch obser- 
vations of K757 by iMeszaros et al] (|2008al ). which show a 
radial velocity shift of 6.2 km s~^, suggest this may be the 
case. If pulsation is present, the "extra absorption compo- 
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Figure 6. Variation of K757 (top panel) and K825 (bottom 
panel) in physical units, as determined from SED fitting. 



nents" visible in the spectra may be pulsation-induced lin e 
doubling, as seen in, e.g., 47 Tucanae (|Lebzelter et al.lboOSi '). 

The location of K757 and K825 on the giant branch, 
along with the variation of these stars in observed colours 
are shown in Fig. [5] and physical units in Fig. (6] It is clear 
from these figures that not only are these two stars among 
the most-evolved in the cluster, but that their variability 
arises from significant changes in both temperature and lu- 
minosity. 

Stellar parameters for K757 and K825 were calculated 
using the same method as for K238 ( t|3.3l). with da t a from 
various optical catalogues, 2MASS and iBover erall \200dt ). 
For K757, this analysis gives a temperature of 4489 ±2017^", 
log(g) = 0.746, and L = 1426 ± 255 Lq for K757; for K825, 
we find r = 4411 ± 155 K, log(g) = 0.662, L = 1615 ± 
227 Lq, whic h arc broa dly consistent with literature values 
(iFrogel et all 119 83: Sn eden et all Il997l : lAlonso et all |200G| : 
IMeszaros et al I f2008bl ). 

Variation in these parameters was estimated by refitting 
the SED using 50-day averages of the optical {g' , V, R, i' , I) 
magnitudes. Fig. [5] shows both stars exhibit a variation be- 
tween a high-temperature, high-luminosity state, and a low- 
temperature, low-luminosity state, with ST^s ~ ±50 K, and 
5L ^ ±30 lJ. (We presume that, since variability declines 
with increasing wavelength, there is negligable variation at 
longer wavelengths than /-band.) If we assume that the to- 
tal optical depth of any dust envelope is largely constant 
through these changes, and that the stellar spectrum is close 
to a blackbody, this corresponds to a radius change of ±1.1% 



Though the absolute errors are large, relative errors can be 
considered much smaller. Ranges also include a small (estimated 
~10 - 20%) noise contribution. 
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around an average of 62 R© for K757, and ±0.9% around 
an average of 69 Rq for K825. 

This type of variation — where maximum brightness 
occurs near-simultaneously for wavelengths longer than ~1 
/xm — is also seen in stronger met al-poor pulsators in 
globular clusters (cf. m C en V42 — iDickens et alj 1 19721 : 
enzies fc Whitelocklll985l V Similar changes can occur on 
spotted stars, though a rotation period of 350 days requires a 
rotation rate of 11 km s~^. This is somewhat large for a star 
this evolved (|Cortes et al.|[2009l ). Furthermore, it would re- 
quire the spots to be relatively constant in size and strength 
over a period of ~1 year; it also cannot explain the changes 
seen in the Ha line. 

A 350-day sinusoidal variation in radius with ~1% am- 
plitude yields a peak pulsation amplitude of only ±0.1 km 
s~^ at the photosphere for both K757 and K825. This value 
is very small compared to both the ~60 km s^^ escape veloc- 
ity, and the speed with which we might expect a dusty wind 
to be accelerated via photon pressure alone. However, the 
speed may increase as the pulsation travels outward through 
the more rarified atmosphere. 

For comparison, we can calculate the speed one might 
expect of a dusty wind. Assuming canonical values of 10 
km s~^ for a 1 000 Lq solar metallicity star, and f cx 
^q0.5[Fo/hi^o.25 dvan LoonI [2OO0I : [Marshall et al.l l2004 ). we 
might expect this speed to be u « 450 ms~^ if gas-dust 
coupling is maintained. If gas-dust coupling is not main- 
tained, and a driving mechanism is coupled to either the 
gas or dust, then we might expect the driven medium to 
escape, and the other to be left behind and fall back onto 
the star. 



4-2.2 K825: Ha line and gaseous mass loss 

Spitzer photometry (|Bover et al.l l200 (j) shows no infrared 
excess around K757 and K825, suggesting a lack of circum- 
stellar dust production. It appears that these stars are of 
too low a metallicity to sustain adequate dust production to 
drive a stellar wind. In the absence of a dust-driven wind, we 
turn our attention to the stellar chromosphere: blue-shifted 
absorption cores, often coupled with line emission wings, 
have been observed in several notable lines and are taken as 
tracers of the mass-loss rate and velocity of an outflow. Such 



ious UV lines fe.E. iDupree et al. 


1984; Duprcc ct al. 1999; 


Cacciari et al. 20041: iDuoree et al 


..2005: Duprce et al...2007.: 


McDonald & van LoonI 20071: Duoree et al.l 2009f). 



We can compare our light curve to recent spectroscopic 
observations of K825 taken as part of an unrelated pro- 
gram me (K757 has already been presented in lMeszaros et ahl 
I2OO9I ). Fig. [7] presents spectra taken with the MIKE double- 
echelle spectrograph mounted on the Magellan/CLAY 6.5 m 
telescope at Las Campanas Observatory using a 0.75" x 
5" slit, giving a resolution oi R ~ 40 000. The spec- 
tra were reduced with bias subtraction, flat-fielding and 
sky subtraction using the MIKE-IDL pipeline updated by 
J. Prochask4f|. Th-Ar arc exposures bracketting the stellar 
targets were used to determine the wavelength scale. Four 



http: / /web. mit.edu/~burles/www/MIKE/ 



epochs are available: two fortuitously near subsequent pho- 
tometric minima, the other two epochs being shortly af- 
ter adjecent photometric maxima. For each spectrum, the 
line bisector is also shown and the velocity shift of the 
line core with respect to the photospheric rest velocity in- 
dicated. These velocities are generally related to the wind 
outflow velocity but do not necessarily approximate them 
modelling suggests that actual out flow velocities are 



higher than indicated by the Hq proflles dMauas et al.|[2006l : 
iMcDonald fc van Loonli2007l : [Mfezaros et al.ll2009t ). 



Variation seen in our lightcurve of K825 is to some ex- 
tent repeated in its Ha and Ca 11 K line profiles: the spectra 
taken at photometric minimum show a moderate asymme- 
try in the Ha emission, while the spectra taken near max- 
imum light show either zero or reversed asymmetry in the 
Ha wings. The Ca 11 K lines show similar variations in the 
strength of their emission, though the red emission com- 
ponent is always dominant here. The behaviour is not en- 
tirely reproducable between epochs — there is considerably 
greater asymmetry at phases 1.32 and 2.16, presumably due 
to a more massive outfiow. This suggests a more chaotic vari- 
ation of the chromosphere, partially but not entirely coupled 
to the pulsation period. 

It is interesting to note that the Ha and Ca 11 line cores 
are permanently blue-shifted. There is no obvious correla- 
tion with pulsation period here either, though the velocities 
of Ha and Ca 11 correlate well with each other. The perma- 
nency of this blueshift, despite substantial changes in the 
emission from the (lower altitude) chromosphere, strongly 
suggests that there is a permanent bulk outfiow from the 
star. As Ca 11 K absorption occurs higher in the atmosphere 
than Ha, and the Ca 11 blueshifts are larger, we can also 
state that the wind is strongly accelerating in this region. 
This means that star's mass loss may be enhanced by pulsa- 
tion, and there may well be other driving mechanisms (e.g. 
magnetic reconnection emitting Alfven waves) which oper- 
ate within the line-forming region of the extended atmo- 
sphere^ 

iMeszaros et al.l (|2009l ') find that K757 shows similarly 
permanent blueshifts in its Ha line cores, despite changes 
in the line emission from the chromosphere. By modelling 
the Ha line, they imply that K757's gaseous mass-loss rate 
changed by a factor of ~6 over their observations. Given 
both stars show a permanent outfiow, with velocities that 
change by a factor of < 2, we may speculate that the wind 
driving mechanism remains relatively constant. The strong 
changes in emission strengths, however, suggest that the rate 
at which material is injected into the wind is strongly vari- 
able. This further suggests that pulsations may vary the 
mass-loss rate from these systems, but that the kinetic en- 
ergy that subsequently accelerates the wind may be provided 
by another, more-constant source. 

4.2.3 K757, K825 and the P ~ L diagram 

As a populous cluster, we may therefore naively expect 
M15 to contain a relatively large number of high-amplitude 
LPVs and SRVs. For comparison, uj Cen — which is five 
times as massive and 0.7 dex more metal rich than M15 — 
has ~20 variables with periods >20 days and F-band am- 
plitudes >0.2 mag llciementlll997l: iMcNamara et al.ll200i : 
iPasauah et~aLl l2004l : Ivan de Ven et al.ll2006l ). By this esti- 
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Figure 7. Variation in the Ha and Ca II K profiles of K825, as observed with Magellan/MIKE. Ha line bisectors are shown as dotted 
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mate, we might expect 4±2 LPVs in M15. Finding two LPVs 
is therefore not unexpected, especiaUy given our incomplete 
coverage. The absence of any high amplitude LPVs in M15 is 
notable, however. This confirms that stellar pulsation plays 
less of a role in stellar evolution at su bstantially sub-solar 
metallicities (|Frogel fc Whitelock|[l998h . but is still present 
among the oldest, most-metal-poor giants we see. 

The long periods of K757 and K825 are unusual. On 
a. P — L d iagram (Fig. El), they fall neat ly onto Wood's se- 
quence D (|Wood et alJI 19991 : IWoodlbOOd ). This sequence de- 
notes pulsators with long, secondary periods (LSPs), though 
we find no shorter-period pulsations that would fit sequence 
C (at ~90 days) or higher harmonics (cf. Fig. [4]). However, 
the two stars also fit sequence E from Derekas et al., which is 
attributed to ellipsoidal variation in red giant (contact) bi- 
naries. Note, however, that K825 shows a more temporally 



asymmetric variability than many similar stars on sequence 
E (cf. 2.5873.24 in Derekas et al.. Fig. 2). For comparison, 
the unconfirmed Bootes I dwarf galaxy member would lie 
near sequence C. 

Pulsation period is general thought to decrease towards 
low metallicity a mon g red giants of similar age in Galactic 
globular clusters (iFrogel fc Whitelock 1998). Were this the 
case, it would be extremely surprising to find K757 and K825 
placed so neatly on one Wood's sequences. Furthermore, 
K825 shows a period as long as the star with the longest 
period in Frogel & Whitelock's sample (NGC 5927 V3, at 
312 days), but at a metallicity wh ere no var i ables this long 
are found; the closest analogue in IClementI l| 19971 ) is NGC 
2419 VIO, a semi-regular variable with P = 81.3 days and 
[Fe/H] = -2.12. Frogel fc Whitelock's data do not include 
all globular cluster variables, and they only claim their cor- 
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relation of pulsation period and metallicity exists for [Fe/H] 
> —1.3. It appears that, while pulsation amplitude declines 
with metallicity, pulsation period does not change apprecia- 
bly. 

In Fig. [9]we reproduce Frogel & Whitelock's Fig. 4 with 
K825 and NGC 2419 VIO included. From these data, it ap- 
pears as though no variables with P > 200 days and [Fe/H] 
< —1 .3 had been discovered prior to iFrogel fc Whitelockl 
l|l998t ). We appear to now have found stars in this miss- 
ing region. The apparent lower luminosity of the RGB tip 
and lack of high-luminosity AGB stars in M15 (M15 has no 
obvious AGB stars ( i)4.2.1ll. in contrast with intermediate- 
meta l licity clusters, e.g ., Ivan Loon et al.ll2006al : I Bover et al.l 
I2OO9I : iMcDonaldl [2009l i means that these stars are on the 
LSP sequence D, rather than appearing in sequence C like 
their metal-intermediate counterparts. 

4.2.4 Intracluster medtum and dust production 

High stellar temperatures and low metallicity make it hard 
to form circumstellar dust around stars in M15. Despite 
this, M15 remains the only cluster i n which intraclus- 
ter gas and dust have been detected |Evans et"aL _ 20031: 
'van Loon et al.1 l2006bl : iBover et al] l200g ~ Ivan Loon et al.l 
2009). These intracluster media must have been produced in 
the (astronomically-)recent past (~1 Myr). However, while 
K757 and K825 both show evidence for gase ous mass loss 
(|Meszaros et al.1 l2008al : iMeszaros et al.ll2009|). ne ither star 
appears to be producing dust (|Boyer et al.l2006l) suggest- 
ing dust production in M15 may not be dominated by stars 
near and above the RGB tip, as seen in other, more-metal- 
rich clusters (e.g. uj Centauri, at [F e/H] ~ -1.6; and 47 Tuc 



I I — 1 

at [Fe/H] = -0.7 — Ivan Loon et a l. 2006a; iLebzelter et al.l 



I2OO6I ; lOriglia et aLll2007l ; iMcDonald et al.ii2009h 



Two likely possibilities remain to explain the presence 
of intra-cluster dust in M15. Du st production in in divid- 
ual stars could b e epis odic (e.g. lOriglia et al. 1 12002| . 2007; 
iMcDonald et al.1 |2009| ). Alternatively, relatively-massive 
(Mh ~0.1 Mq) amounts of dust-containing medium could 
be produced in a superwind at the AGB (and possibly 
RGB) tips, or in creating a post-AGB star. Indeed, Boyer et 
al. observe dust in M15's planetary nebula, Pease 1. How- 
ever, the amount of in tra-cluster medium (Mh ~0.3 Mq; 
Ivan Loon et all l2006bl ) would require this short phase to 
have recently occurred near-simultaneously (within ~1 Myr) 
for at least three stars. As a result, we expect that episodic 
dust production around stars with more-sustained gaseous 
outflows seems the more-likely explanation for the presence 
of intra-cluster dust, though both factors could contribute. 

If dust production is episodic, would we expect to see 
dust in the most-evolved stars? Of the 17 most-evolved stars 
in M15 (Fig. [5]), only K479, K421 a nd K373 show clear in- 
frared excess in iBover et al] l|2006l ). K204 may also show 
excess, but it is blended with K570 in the Spitzer Space 
Telescope 3-70-/im photometry. Assuming these stars have 
luminosities of 1400-1 600 Lp) (based on ^4.1.711 . we can com- 
pare this to Fig. 20 of IMcDonald et aU l|200g ). which shows 
the fraction of dusty stars as a function of luminosity in uj 
Cen. Assuming a detection limit for circumstellar dust cor- 
responding approximately to the 6cr line in McDonald et al. 
(this approximately takes into account the metallicity dif- 
ference between the clusters), we might presume that ~30% 
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(five) of these 17 stars would have circumstellar dust, which 
is statistically indifferent from the three or four observed. 

The lack of dust production in K757 and K825 is puz- 
zling, however, when one considers th at the above Ha obser- 
vations ( ^4:2:21 iMeszaros et al.ll2009^ indicate gaseous mass 
loss does occur from these st ars, and indeed from stars fur- 
ther down the giant branch l|McDonald fc vaii Loonll2007l : 
iMeszaros et alj 20091 ) . The pulsation period is long enough 
for dust to form between pulsations, but not long enough 
for dust to become so cool as to be undetectable in its [8]- 
[24] colour (assuming Vac ~ 50 km s~^). The pulsations also 
appear not to generate stron g shock waves (which show u p 
as sharp peaks in Ha — e.g. iMcDonald fc van Loonll2007l ). 
which may assist dust production. The constant, moderate- 
velocity outflow provided by the chromosphere may also pre- 
vent the formation of a quasi-static molecular layer required 
to seed dust formation (TsuU et al. 1997); the formation of 
such a layer being additionally hindered by the hotter pho- 
tospheric temperatures and reduced availability of metals in 
such low-metallicity stars. Furthermore, it may be that the 
low amplitude of the pulsations means the scaleheight of the 
atmosphere is not sufficiently raised, meaning the wind den- 
sity beyond the dust-formation radius is too low to efficiently 
produce dust, while the increased effective temperature of 
the metal-poor stars means the dust-formation radius is too 
far from the star. 

It i s clear from observation s of more metal-rich 
clusters llRamdani fc J orisscn 2 00ll: Ivan Loon et al. I l2006al : 
iLebzelter et al.1 l2006l : lMcDonald et al.1 120091 ') that eflicient 
dust production occurs mainly in strongly variable stars at 
the upper end of the RGB, and that the amount of mass 
loss roughly correlates with the strength of the pulsation 
as traced by the optical photometric amplitude. Our obser- 
vations here confirm that correlation, finding no detectable 
dust production in stars which (while at the RGB tip) ex- 
hibit only weak pulsation. 



5 CONCLUSIONS 

This study presents a search for long-period variability 
among giants in the globular cluster M15 and has confirmed 
the most metal-poor variables known in our Galaxy: K757 
and K825. These stars are very close to the RGB-tip and 
show lightcurves characteristic of pulsation, albeit with very 
low amplitude. Their periods pla ce them on the 'long sec- 
onda ry period' (LSP) sequence l|Wood et al.l 1 19991 : IWoodI 
|2000D . though no 'primary' period has been found. 

The pulsational velocities are similar to those measured 
at the base of the wind, but it is not certain that pulsation is 
required to launch or drive the wind. In any case, these pul- 
sations appear not to lead to dust production, which must 
be caused by a different mechanism, possibly episodically. 
Although we do find evidence for acceleration in the singly- 
ionised calcium line profiles, radiation pressure on grains 
cannot therefore be held responsible for driving the wind. 
Instread, Alfven waves may couple to a weakly-ionised gas 
and thus drive the wind, something which cou ld be facili- 
tated at lower metallicity (jvan Loon et al.ll2010l ). 
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